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Nonuniform background electromagnetic fields, once implemented in lattice quantum chro¬ 
modynamics calculations of hadronic systems, provide a means to constrain a large class of 
electromagnetic properties of hadrons and nuclei, from their higher electromagnetic moments 
and charge radii to their electromagnetic form factors. We show how nonuniform fields can 
be constructed on a periodic hypercubic lattice under certain conditions and determine the 
precise form of the background 17(1) gauge links that must be imposed on the quantum chro¬ 
modynamics gauge-field configurations to maintain periodicity. Once supplemented by a set 
of quantization conditions on the background-field parameters, this construction guarantees 
that no nonuniformity occurs in the hadronic correlation functions across the boundary of 
the lattice. The special cases of uniform electric and magnetic fields, a nonuniform electric 
field that varies linearly in one spatial coordinate (relevant to the determination of quadru¬ 
ple moment and charge radii), nonuniform electric and magnetic fields with given temporal 
and spatial dependences (relevant to the determination of nucleon spin polarizabilities) and 
plane-wave electromagnetic fields (relevant to the determination of electromagnetic form 
factors) are discussed explicitly. 


I. INTRODUCTION 

The electromagnetic (EM) properties of hadrons and nuclei arise from the underlying interactions 
among their strongly interacting quark constituents and the photons. Consequently, any reliable 
theoretical approach that aims to elucidate the electromagnetic structure of hadronic systems must 
necessarily account for the nonperturbative nature of quantum chromodynamics (QCD). Lattice QCD, 
which is a Monte Carlo evaluation of the QCD path integral regulated through a hnite discrete space- 
time, is the only such method by which to perform hrst-principles calculations of hadronic systems. It 
relies on the fundamental degrees of freedom of QCD, i.e., quarks and gluons, and as input, only takes 
the QCD parameters, the mass of quarks and the strength of the coupling constant. To be compared 
with nature, the results of lattice QCD calculations must be systematically extrapolated to the contin¬ 
uum and inhnite-volume limits. Quantum electrodynamics (QED) can be treated along with QCD in 
a similar fashion, although the computational cost associated with a Monte Carlo evaluation of QCD 
and QED path integral is considerably higher, see Refs. mu for recent progress in this direction. 
Alternatively, given that EM interactions are perturbatively small, one can constrain the EM structure 
of hadrons and nuclei, for example their EM form factors, by evaluating the matrix elements of current 
operators in the presence of solely QCD interactions. This approach has met with success in several 
cases, but is generally challenging (see Ref. m for a review of lattice QCD calculations of nucleon 
structure). It is therefore useful to consider further alternatives. 

When interested in the response of hadronic systems to weak EM helds, which determines properties 
such as EM moments and polarizabilities, a powerful method is to introduce electromagnetism through 
a classical background held. In this approach, hxed U(l) gauge links are simply imposed on the 
QCD gauge links in the lattice formulation. Depending on the computational resources available and 
the type of quantities that are being considered, this imposition may be performed solely on the 
valence quark sector of QCD - where only the computation of quark propagators is inhuenced by the 
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additional U{1) links - or on both the sea and valence quark sectors - where the U(l) gauge links 
are also incorporated in the generation of the QCD gauge configurations. The former is not a reliable 
approximation when there are sea-quark (disconnected) contributions to hadronic correlation functions 
in background fields. A low-energy (multi-)hadronic theory that describes the interaction of the hadron 
(or nucleus) with weak external EM fields can be matched onto appropriate lattice QCD correlation 
functions. This matching constrains those parameters of the hadronic theory that characterize the 
response of the hadron (nucleus) to the applied external fields |13H24| . While uniform background 
fields already provide a means to constrain a wealth of EM quantities such as the magnetic moment 
and the electric and magnetic polarizabilities (see Refs. [22l|23l|25H37] for recent progress in discerning 
EM properties of hadrons and light nuclei using lattice QCD with the background field method), more 
general background fields potentially provide sensitivity to many additional parameters of the low- 
energy theory. In particular, spin-dependent structure parameters beyond the magnetic moment can 
show up at low orders in the weak-field strength if one allows the spatial and/or temporal derivatives of 
the background fields to be nonvanishing. Among such quantities are spin polarizabilities of nucleons 
|221l38j . and the quadrupole moment of hadrons and nuclei with spin, s > ina. Moreover, nonuniform 
background fields do not constrain one to the static limit of EM form factors, and by injecting energy 
and momentum into the system, provide a means to directly evaluate the corresponding off-forward 
hadronic matrix elements of the EM current from a response to the background fields |39| , with different 
systematic uncertainties than other methods. Another application of such nonuniform background 
fields, as is recently proposed and implemented in Ref. 00], is to evaluate the hadronic vacuum 
polarization function (as the leading hadronic contribution to the muon anomalous magnetic moment) 
through evaluating the magnetic susceptibilities with a plane-wave background field. 

In order to explore such possibilities, one first needs to properly implement the desired EM back¬ 
ground fields in lattice QCD calculations. A class of nonuniform EM fields, for example, have been 
implemented in Refs. 02021 to obtain some preliminary results for spin polarizabilities of the nucleon 
|42) . These studies do not consider background gauge potentials that are periodic at the boundary. 
Retaining the periodic boundary conditions (PBCs) that are imposed on the gauge fields in the major¬ 
ity of lattice QCD calculations requires a nontrivial implementation of the U{1) gauge links. However, 
such an implementation enforces smooth behavior of correlation functions across the boundary of the 
lattice and is a desired feature. Moreover, quantifying the behavior of the system in the hadronic 
theory is generally more straightforward in periodic EM potentials. Although nonperiodic implemen¬ 
tations of uniform background fields have been pursued in some earlier lattice QCD studies (by placing 
hadronic sources away from the boundary effects) |25H30) . quantifying uncertainties associated with 
these nonuniformities is difficult 123031221 . These issues can be prevented by explicitly modifying the 
naive 17(1) gauge links, and imposing conditions on the background field parameters, such that when 
setting the value of the gauge links at one boundary of the lattice to its value on the opposite bound¬ 
ary, no nonuniformities occur in the value of the 17(1) plaquettes. This guarantees smooth behavior 
of hadronic correlation functions across the boundaries. This paper presents in detail a procedure for 
the periodic implementation of background gauge fields, under certain conditions that are enumerated, 
along with several examples that follow from these general considerations. 

A condition that supports periodic background 17(1) gauge fields on a hypercubic lattice^ is well 
known for the case of uniform EM fields, namely the’t Hooft quantization condition (QC) |13[ ITTl I25T 
07 ], and has been commonly implemented in lattice QCD calculations with the use of background 
fields |22[ l23l l3Tll37[ 148] . This condition requires the magnitude of the electric, E, or magnetic, B, 
field on a torus to be quantized, and follows from a simple argument: for a closed surface geometry. 


^ Although the words “hypercubic” and “hypertorus” are used throughout to refer to the lattice geometry, this paper 
considers the more general case of an anisotropic geometry where both the lattice spacing and the extent of the volume 
in temporal and spatial directions are different. In particular, the general result in Sec. |III| accounts for distinct lattice 
spacings and volume extents in all directions. 
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the net flux of the EM field through that surface is required to be quantized. The same condition can 
be obtained by imposing a more general boundary conditions, namely electro/magneto-PBCs. These 
boundary conditions require the gauge and matter fields to be periodic up to a gauge transformation, 
and have been studied by ’t Hooft for the case of Abelian and non-Abelian gauge theories, to explore the 
properties of the flux of the corresponding field strength tensors in the confinement regime |45l SHI |50] . 
When spacetime is discretized on a hypertorus, special care must be given to the implementation of 
the 17(1) gauge links at each spacetime points, see e.g.. Refs. |48l ISTl 152) . to guarantee the values of 
the 17(1) elementary plaquettes remain constant in uniform background fields. 

To implement background EM fields with arbitrary space-time dependences, a similar procedure 
must be undertaken to ensure that, despite the gauge links having been set to satisfy the PBCs at the 
boundary, the expected values of the 17(1) elementary plaquettes are correctly produced adjacent to the 
boundaries. This requirement enables one to determine the modified links near the boundary, as well as 
the conditions that the parameters of the chosen background fields must satisfy. We first demonstrate 
this procedure for a rather special case of an electric field generated from a scalar gauge potential 
with an arbitrary dependence on only one spatial coordinate. This allows one to generalize, rather 
straightforwardly, to the case of background fields generated from more general gauge potentials, but 
also makes one appreciate the subtleties and limitations encountered in the general case. We use this 
special case to study in detail the examples of a uniform electric field and an electric field that varies 
linearly in one spatial coordinate. We verify the periodicity of the setup by numerically evaluating 
the correlation functions of neutral pions in these background fields. This special case, along with the 
examples, is presented in Section [TT| The general considerations are presented in Section m where 
it is shown that a periodic implementation of gauge links in our framework is possible if the flux of 
the electric and magnetic fields through each plane on the lattice is independent of the coordinates 
transverse to that plane, allowing the flux to be quantized. The same condition also arises from 
consideration of the functions introduced in the modified links adjacent to the boundary such that 
the expected value of plaquettes are produced in all planes on the lattice. We take advantage of our 
general results, as summarized in Eqs. (51 53) of this paper, to work out several phenomenologically 
interesting examples. These examples focus on the background fields that give access to some of the 
spin polarizabilities of nucleons as suggested in Ref. |22| . as well as the case of a plane-wave EM fields 
as proposed in Refs. |39| 140) . These examples are presented in Section IV We summarize our results 
and conclude in Section |V] An appendix is devoted to demonstrating the connection between our 
results concerning background 17(1) gauge fields with PBCs and those obtained under the imposition 
of certain electro/magneto-PBCs. 


II. A SPECIAL CASE: AN ELECTRIC FIELD ARISING FROM GAUGE POTENTIAL 

= (Ao(x3),0) 


Let us choose a periodic 17(1) gauge field 


— (Aq, —A) — ( Ao(x3 — 


1 L J 


L ),0 


( 1 ) 


with an arbitrary dependence on the third component of the position three-vector, X 3 , such that a 
periodic electric field is generated in the X 3 direction. 


E — —V Aq — E/(x3 


L 1 


L)x3. 


( 2 ) 


With the use of the floor function in the argument of the functions, the fields are made periodic in X 3 
with periodicity L. As is seen, we have adopted a mostly negative signature for the Minkowski metric. 
Throughout this paper, we take any boldfaced letter to denote a three-vector. We leave the letters 
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boldfaced even if they correspond to the components of a three-vector as to distinguish them from the 
components of a Minkowski four-vector, e.g., Xj = —Xj for = (t, —x). 

The background field is implemented in a lattice QCD calculation by multiplying the QCD gauge 
links by the U{1) gauge links through a direct product. Explicitly, for the choice in Eq. Q, 




^ gieQAo(x3-["^]L)atX<5^,0 


(3) 


where Q denotes the electric charge operator. Here, and in what follows, we define and at to denote 
the lattice spacings along the spatial and temporal directions of the lattice, respectively. For a periodic 
lattice with spatial extent L and temporal extent T, the value of the plaquette with 0 < X 3 < L — 
and 0 <t <T — at in. the 0 — 3 plane is 


^(0,3)(x3,t) = 
X 


f^WCD)t(x 3 ^ 


— g*eQMo(x3)-Ao(x3+as)]at (X 3 t) 


( 0 , 3 ) '.^3,'^;, (4) 

where we have left implicit the dependence on the xi and X 2 coordinates. Note that in the continuum. 


d> 


(E) 

( 0 , 3 ) 


(X 3 ) = 



E3(x3)dx3 


[Ho(x 3 ) - Ho(x 3 -I- Os)] at, 


(5) 


is nothing but the total electric flux through the surface area of the elementary plaquette originated 
from the point (x 3 ,t) in the 0 — 3 plane. Therefore, the desired value of the plaquette is obtained in 
Eq. (Q. 

Since the lattice action depends upon links that originate from points Xj = L — lort = r—1 and 
end at points Xi = L or t = T, one is required to specify the boundary conditions. By choosing PBCs, 
we demand that the f7(l) gauge link to be periodic according to Eq. (Q. Then the value of the link 
that originates from X 3 = L is set equal to its value at the origin. As a result, one must examine more 
carefully the value of the plaquettes located at X 3 = L — a^. 


P(o,3)(i-a.,t) = U^'^^^\L-as,t) t/f - a,, t + a*) 

X (0,t) - as,t) 

— ^ieQ[Ao{L-aa)—Ao{L)jat ^ieQ[Ao{L)-Ao{0)jat 


( 6 ) 


where Aq is defined as the de-periodified counterpart of Aq, Aq = Ao((x 3 ~ [x]) + [¥] )■ Then the first 
phase factor correctly accounts for the total electric flux through a plaquette located at X 3 = L — Ug, 
while the last phase factor must be eliminated. To achieve this, we are free to introduce an additional 
link in the X 3 direction. 


^(QCD)^^^ ^(QCD)^^^ ^ g*eQatAo(x3 —[^]l,)x5^_o ^ g*eQ[j4o(0)—ylo(i)](t—[;f;]T)X(5^_3(5x3,i- 


(7) 


such that the last phase is canceled from Eq. ([^ . Clearly, the additional link does not affect the value 
of the adjacent plaquettes. 

Since the additional link introduced in Eq. ([^ is t dependent, and given that the gauge link is also 
required to be periodic with respect to the time variable, one must study the value of the plaquette 
located at X 3 = L — and t = T — at more closely. For this plaquette, which is located at the far 
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corner of the lattice,^ 

P(o,3)(i - a,,T-at) = - at) - a„ 0) 

X at) - as,T - at) e-*'=‘3[^o(o)-Ao(L)](T-at) 

^ ^ieQ[Ao{L-as)-Ao(L)]at ^(QCD)^^ - as,T - at) X (8) 

Although the first phase factor correctly accounts for the total flux through a plaquette located at 
X 3 = L — as, the last factor modifies the value of this plaquette away from the desired value. By 
demanding 


^ieQ\A(,{0)—Ao{L)]T _ ^ 


(9) 


the value of this plaquette is fixed to its desired value. Eq. ^ places a QC on the net flux of the 
electric field through the 0 — 3 plane on the lattice, 




(E),net _ 

( 0 , 3 ) = 


Ao( 0 )-Ao(L) T = eQ 


dt 


E3(x3)dx3 = 27 rn, 


( 10 ) 


where n G Z. This QC could indeed be deduced a priori by recalling that the 0 — 3 plane represents 
a closed surface area due to PBCs (the surface area of a torus), and the net flux through this surface 
must necessarily be quantized. Eqs. ([^ and (10) are the main results of this section. Two particular 
cases of a uniform electric field and a linearly varying electric field along the X 3 direction can be readily 
worked out from this general result, and are presented in the following. 


Example I: A constant electric field in the X3 direction 

An external periodic U{1) gauge field 


Afj,= (-E X {x3- R- 


L J 


^),o). 


with constant E and R, gives rise to a uniform electric field in the X 3 direction, 

E = EX3. 

This electric field can be achieved by implementing the link, 

jy(QCD)^^^ ^(QCD)^^^ ^ g—ieQEat(x3 —R—[^]L)x5^,o gjeQEL{t—[|,]T)X(5^_3(5x3,i-a 

which gives rise to the desired value of the plaquette in the 0 — 3 plane 

P(0,3)(x3,t) = V'§^^\x,,t), 


( 11 ) 


( 12 ) 


(13) 


(14) 


for 0 < X 3 < L — Os and 0 < t < T — at- In addition to the modification to the naive 17(1) links as 
presented in Eq. (13), the plaquettes at the boundary are correctly valued if the magnitude of the 

E=^A^, (15) 


electric field is quantized. 


eQTL 


^ Strictly speaking, there is no notion of a corner in a toroidal geometry. However, once an origin is specified for the 
coordinate system, we can choose to define the corner as the plaquette located at X 3 = 1/ — as and t = T — at- 
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FIG. 1: A comparison of choices of background gauge fields that result in a uniform electric field in the 
X 3 direction: “Modified links - Quantized” denotes the choice of the gauge links in Eq. (13) with i? = 0, 
supplemented by the QC in Eq. (151, with n = 3. We have quantized the field for the d quark, corresponding to 
IQI = which guarantees the quantization of the field for the u quark as well. “Modified links - Nonquantized” 
corresponds to the same choice of the gauge links but with a nonqunatized value of the field, n = e. “No 
modified links - Quantized” and “No modified links - Nonquantized” denote the naive choice of the gauge links 
without including the additional links in the X 3 direction, with quantized and nonquantized values, n = 3 and 
n = e, respectively. The upper panel shows the correlation function (projected to zero transverse momentum), 
C(x 3 ,t), as a function of X 3 — at a fixed Euclidean time, r/at = 18, while the lower panel depicts the 

dependence of the quantity M(x 3 ,r) = log ~ at r/at = 18. In this demonstrative study, 

the lattice volume is taken to be E = (12as)^ x (24at) with as = at ^ 0.145 [fm], the QCD gauge configurations 
are quenched, the C/(l) gauge links are solely imposed on the QCD gauge links in the valence quark sector, and 
only the connected pieces of the correlation functions are evaluated for the neutral pion. The values of X 3 and 
Xg*^'^'^^ in the figure are in units of a^. The dashed lines denote the boundary of the lattice. 


with n G Z, in agreement with previous works, e.g., Refs. [isKiailelllTlEIllsg. 

We have numerically verified that, in contrast with the nonperiodic case where the additional link in 
the X 3 direction is not implemented, the choice of the background gauge field in Eq. (13), supplemented 
by the QC in Eq. ( |15[ ), results in the periodicity of the correlation functions of neutral pions. The 
upper panel of Fig. llTshows the pion correlation function (projected to zero transverse momentum), 
C(x 3 ,r), as a function of X 3 — Xg at a fixed Euclidean time, r/at = 18 (r = it). Explicitly, this 
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correlation function is defined as 


= C'(x3,r;x^®''‘'\0) 


P1=P2=0 


(-Lg 

Y. E C(x,t;x(”),0), 

xi=0 X2=0 


(16) 


where 


C(x,r;x(-‘^),0) = (0|a(x, r)Ot (x(-‘^), 0)|0 )e • (17) 

oIIOtt) is a lattice interpolating operator that creates (annihilates) any hadronic states with the 
quantum numbers of the neutral pion. Subscript E refers to the fact that the expectation value is 
evaluated in the background of an electric field, E. The calculation only involves imposing the 17(1} 
gauge links on the QCD gauge links in the valence sector, denotes the location of the source, 

which for the upper panel is taken to be = (0,0,0). Since for a neutral pion in a uniform 

electric field, the finite-volume correlation function with PBC must be symmetric about the point 
j -j- the deviation of the correlation function from symmetricity for nonperiodic gauge-link 

choices, including those with the correct link structure but with nonquantized values of electric field, 
signals the breakdown of translational invariance in units of L in the X 3 direction (this translational 
invariance is the analogue of the magnetic translation group discussed in Ref. m for a uniform 
magnetic field). Such breakdown is most evident in the quantity 


M(x 3 ,t) = log 


C(x3,t) 
C'(x3,r-t- 1)’ 


(18) 


as is plotted as a function of X 3 — for r/a^ = 18 in the lower panel of Fig. Here, the source is 
located = (0,0, Da^) and therefore the boundary point X 3 = L = 0 corresponds to X 3 — 

in these plots. Nonuniformities in M(x 3 , r) when crossing this boundary (denoted by the dashed line) 
are observed in all the cases considered, except for the “Modified links - Quantized” case.^ This is 
again a signature of losing translational invariance in units of L in the X 3 direction. In Refs. |231 HH] . 
a similar kinked feature was observed in the correlation function of neutral pions with nonperiodic 
implementations of a uniform electric field with the choice of a time-dependent gauge potential. 


Example II: A linearly varying electric field in the X3 direction 

An external periodic U(l) gauge field 

En 


A,= [-^(x3-R-[-^\ L)^0 
gives rise to a linearly varying electric field in the X 3 direction, 

E = Eox(x3-R- ^ L)x3, 


(19) 


( 20 ) 


as plotted in Fig. This electric field can be implemented in a lattice QCD calculation through the 
following links. 




® These nonuniformities may be quantified more precisely by evaluating (the finite-difference approximation to) the 
derivative of the functions with respect to X3. As the continuum limit is approached, this (numerical) derivative 
diverges near the boundary as a result of nonperiodic implementations. 










FIG. 2: The scalar potential in Eq. ( 19 ) (the left panel) with Eq < 0 and i? = ^ is a finite harmonic oscillator 
potential between 0 < X 3 < L. It produces a linearly varying electric field, Eq. (20), as depicted in the right 
panel. The periodic images of the potential and the electric field are also shown in the figures. 


which gives rise to the following value for a plaquette in the 0 — 3 plane 




(0 


,3)(x3,t) = e-Q^oa.a.(x3-it+^) p(QCD)^ ^ 


( 22 ) 


for 0 < X 3 < L — Os and 0 < t < T — at- There exists a QC when R ^ j which constrains the value of 
the slope of the electric field to be quantized, 


En = 


2Trn 


eQTL{-R+^) 


(23) 


with n £ Tj. When the offset value R is chosen to be , Eq. (^9^ trivially holds and no quantization 


constraint is placed on Eq.'^ Moreover, the additional link at point X 3 = L — in Eq. (21) will be 
equal to unity. 

We have implemented a linearly varying background electric field, with and without the additional 
links in the X 3 direction, see Eq. (21), for quantized and nonqunatized values of the electric field 


slope, see Eq. (23). Fig. [^corresponds to the choices of electric field with the offset value R = 0. 
Consequently, according to Eq. (23), a QC is necessary to guarantee the full periodicity. In the left 


panel of the figure, the neutral pion correlation function (projected to zero transverse momentum), 
C(x 3 , r), is plotted as a function of X 3 — at a fixed Euclidean time, t jat = 18. For a neutral pion 
in an electric field that varies as a function of xg, the finite-volume correlation function with PBCs 
will no longer be symmetric about the point ^ -|- Xg'^'^'^^ (see Ref. jl5) and discussions associated with 
Fig. 0 - a feature that is indeed observed in the left panel of the figure. However, by displacing the 
source to Xg^'^'^^ 7^ 0, the breakdown of translational invariance in units of L is again manifest through 
nonuniformities when crossing the boundary of the lattice. This feature, as is plotted in the right panel 
of the figure, is observed in quantity M(x3,r) at a fixed time when crossing xg — Xg^’^'^^ = 30 ^, with 


When there is no QC placed on the slope of the field, Eq, the field can become arbitrarily strong at the boundaries of 
the lattice in the X 3 direction in the large-volume limit. When interested in the response of the system to the external 
field at leading orders in the held strength, one needs to make sure that the held remains weak at any spacetime point 
on the lattice. This can be achieved by tuning Eo to have appropriately small values, e.g., values that are power-law 
suppressed in volume. 
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FIG. 3: A comparison of choices of background gauge fields that result in a linearly varying electric field in 
the X 3 direction: “Modified links - Quantized” denotes the choice of the gauge links in Eq. (21) with R = 0, 


supplemented by the QC in Eq. (23), with n = 3 and \Q\ = “Modified links - Nonquantized” corresponds to 
the same choice of the gauge linksBut with a nonqunatized value of the field slope, n = e. “No modified links - 
Quantized” and “No modified links - Nonquantized” denote the naive choice of the gauge links without including 
the additional links in the X 3 direction, with quantized and nonqunatized values of the field slope, n = 3 and 
n = e, respectively. The left panel shows the correlation function (projected to zero transverse momentum), 
C(x 3 ,t), as a function of X 3 — at a fixed Euclidean time, r/at = 18, while the left panel depicts the 

dependence of the quantity M(x 3 ,t) = log on X 3 — at r/a* = 18. The details of this numerical 

study are the same as in Fig. The values of X 3 and Xg^'^'^^ in the figure are in units of Ug- The dashed lines 
denote the boundary of the lattice. 


x^src) _ more prominent for implementations that do not use the modified values of the 

gauge links near the boundary. 

According to Eq. (21), when the field offset value is set to i? = j, no additional link is needed 
to guarantee periodicity. As a result, our implementations of both cases are trivially identical as is 
shown in Fig. However, a nontrivial check is to show that the periodicity of correlation functions is 
retained with arbitrary (nonqunatized) values of the electric field slope. Although there exists no QC 
for this case, we have continued to label the choices of the electric field slope as in the case of R = 0. 
Then, as is evident from the right panel of Fig. the function M(x 3 ,r) smoothly crosses over the 
boundary of the lattice (up to the uncertainty of the data points), and correlation functions respect the 
periodicity of the fields. Additionally, the left panel of the figure suggests that the correlation function 
is approximately symmetric about the point ^ + X 3 when R= in contrast with the case of i? = 0. 
This feature is understood by recalling that for a neutral pion, the leading dependence of the correlation 
function on a nonuniform external field arises from its charge radius. The corresponding contribution 
then scales as ~ Eq for the case of a linearly varying field in Eq. (12), which is constant. The next 
contribution results from the nonvanishing polarizability of the pion and scales as E^. Although this 
contribution is X 3 dependent, it only depends on the square of the field, a quantity which is symmetric 


about X 3 = ^ if the field offset value is set to R = ^ (with = 0), see Fig. |2 This also explains 


rn 




the asymmetry in the neutral pion correlation function in the background of a linearly varying field 
with R = 0, for which the field squared is not a symmetric function of X 3 around the midpoint of the 
lattice in the X 3 direction. 
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FIG. 4: A comparison of choices of background gauge fields that result in a linearly varying electric field in 
the xg direction, as described in the caption of Fig. but with the offset value being set to i? = The 
left panel shows the correlation function (projected to zero transverse momentum), C'(x 3 ,t), as a function of 
xg — at a fixed Euclidean time, t jat = 18, while the left panel depicts the dependence of the quantity 
M(x 3 ,r) = log ™ ^3 ~ at r/ot = 18. The details of this numerical study are the same as in 

Fig. 0 The values of xg and Xg^'^'^^ in the figure are in units of Og. The dashed lines denote the boundary of the 
lattice. 


III. A GENERAL CONSIDERATION: ELECTRIC AND MAGNETIC FIELDS ARISING 
FROM GAUGE POTENTIAL = (Ao(x,t),-A(x,t)) 

Having gained experience with the rather special case of an electric field generated from a scalar 
potential that only depended on a single space coordinate, we are ready to consider the case of a gauge 
potential with arbitrary dependences on all space and time coordinates. As will be demonstrated 
shortly, a periodic implementation of 17(1) gauge links can be achieved through modifying links adjacent 
to the boundary for a class of background fields that do not generate coordinate-dependent flux through 
any spacetime plane on the lattice (a closed surface on a torus). 

In order to generate electric or magnetic fields with arbitrary dependences on space and time coor¬ 
dinates, one needs to consider a generic gauge potential with all four components being nonvanishing 
and having arbitrary space-time dependences, 

A^(x,t) = (Ao(x,t),-A(x,t)). (24) 

All the functions are to be understood to be periodified according to the floor-function prescription 
introduced in the previous section. Explicitly, all the Xj dependences of the functions must be under¬ 
stood as a dependence on (x* — L). Similarly, all the time dependences are to be replaced with 

(t — [|;] T) dependences. As introduced above, once we de-periodify the functions in either space or 
time coordinates, we place a tilde over them. If 0 < Xj < L and/or 0 < t < T, there is no differ¬ 
ence between functions with and without the tilde, and in the following this distinction only makes a 
difference when the functions are to be evaluated at Xj = L and/or t = T. 

In the covariant formulation, the electric and magnetic fields are closely related and one does not 
need to distinguish the 0 and i components of vectors (where i = 1,2,3 refers to spatial indices). 
However, to carefully account for different lattice spacings and volume extents in temporal and spatial 
directions, the plaquettes in the 0 — 7 and i — j planes, and consequently the cases of electric and 
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magnetic fields, are studied separately. The goal is to generate the following electric and magnetic 
fields 

E(x,t) = -^^^-Vylo(x,t), 

B(x, t) = V X A(x, t), (25) 


in the continuum limit by imposing the 17(1) gauge links on the QCD gauge links. Naively, one can 
start with the following links 


f^WCD)(x^^) ^ ^ gieQyl,(x,t)a.^ 


(26) 


which then must be supplemented by additional 17(1) phase factors adjacent to the boundaries of the 
lattice to ensure periodicity of gauge-invariant quantities. To find these additional phase factors, it 
suffices to study the value of elementary plaquettes throughout the lattice. We first note that the 
naive implementation of the 17(1) gauge links gives the correct value of plaquettes everywhere except 
at Xi = L — ttg and t = T — at- Explicitly, for 0 < Xj < L — a<j and 0 < t < T — at, the value of 
plaquette in the 0 — i plane reads 


_ ^-ieQ[Ai{xi,t)-Ai{xi,t+at)]as+ieQ[Ao{xi,t)-Ao{xi+as,t)]at ^ t) 


(27) 


where we have dropped the Xj and x^ dependences of the functions for brevity as they are fixed in this 
expression. Here, and in the following, i, j and k assume distinct values. The 17(1) phase appearing 
in the value of the plaquette is clearly the expected value by noting that the electric flux through the 
corresponding surface area in the 0 — i plane is 


rXi-\-as rt-\-at 

= / dx't dt''Ei{x'i,Xj,Xk,t') 

’ J Ki Jt 

= - [Ai{xi,xj,xk,t) - Ai{xi,xj,xk,t + at)] a^ [Ho(xj, x^, x^, f) - Ao(xj -k as,Xj,Xk,t)] at 

-\-0 (^a^at, aftts) ■ (28) 


I'-pN 

In the continuum limit, the flux per unit area, ]^/asat, is finite and nonvanishing for any nonvanish¬ 


ing value of the electric field. As a result, the exponent of the 17(1) phase factor in Eq. (27) correctly 
accounts for the flux of electric field through the corresponding surface area in the continuum limit. 
Since the value of the gauge links at the boundary will be set to their value at the origin due to the 


use of periodic functions in Eq. (24), the value of the plaquettes adjacent to the boundary of the lattice 


will not necessarily produce their value as expected for the electric-field flux through these plaquettes. 
This can be easily seen by evaluating the value of plaquettes at Xj = L — a* and 0 < t < T — at, 

iP(o,i)(xi = L-as,t) = 

X t)e-ieQAo(0,t)at ^y-ieQM{L-aa)as 

— ^-ieQ[Ai{L-as,t)-Ai{L-as,t+at)\as+ieQ\Ao{L-as,t)-Ao{L,t)]at x — a t) 

^^ieQ[Ao{L,t)-Ao{0,t)]at 


where again the Xj and x^ dependences of the functions are suppressed. One may now choose to modify 
the value of the link along either the 0 or the i directions in such a way to cancel out the extra phase 
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in Eq. (29). However, it is easy to show that by modifying the link along the 0 direction, the values of 


all the adjacent plaquettes are affected. In particular, one can not consistently make this modification 
throughout the lattice without shifting the value of another set of plaquettes from their desired values. 
However, one can make the modification along the i direction without changing the value of adjacent 
plaquettes. This is exactly the case we considered in the previous section, Eq. ([^ (upon setting Ai = 0 
in the current case). However, in contrast to the case considered previously, the scalar potential that 
appears in Eq. ( |29[ ) carries an arbitrary t dependence. As a result, the prescription of Eq. ([^ will not 
work and we need to find a more general modified link in the i direction to cancel the last phase factor 


in Eq. (29). It is easy to see that the following modification to the Ui link at x* = L — 


[/(QCD)^^ ^(QCD)^^ ^ ^ieQAi{yi,t)as ^ gie(3[Ao(xi=0,Xj,Xfc,t)-Ao(xi=L,Xj,Xfc,t)]/i,o(xj,xj,,t)X(5x;,L-a 


(30) 


achieves the desired result provided that the boundary function /j^o(xj) t) introduced above satisfies 
the following relation 

Ao(xi = 0,Xj,Xfc,t + at) - Ao(xi = L,Xj,Xfc,t + at) /i,o(xj,x^, t + at) 

Ao(xi = 0,Xj,Xfc,t) - Ao(xj = L,Xj,Xfc,t) (/i,o(xj,x^, t) + a*). (31) 

Note that Xj and x^ denote coordinates that are transverse to the 0 — z plane. In order for the function 
fi^o to not spoil the periodicity of gauge links in the temporal and spatial directions, its dependences 


upon space and time coordinates must be understood through the floor-function prescription. Eq. (31) 


is a recursive relation, and once the initial value of the function, corresponding to /i,o(xj,Xfc,t = 0), is 
input, its value at every other point t, fifi{xj,Xk,t), can be obtained. For the special case of the previous 
section, ft^ is only a function of t and satisfies a simple recursive relation, fifi{t+at) = fifl{t)+at. Once 
fi,o{0) is set to zero, the solution to this equation simply is fifl{t) = t, as already prescribed in Eq. ([^. 
In general, if the function /j^o depends on xj and x^ for the chosen gauge field, additional conditions 
must be placed on /j^o- We will deduce these relations once we extend the above considerations to the 
plaquettes in other planes. 

The fifi function, as well as other gauge fields, have explicit time dependence and are periodic in 
this variable. As a result, fixing the value of plaquettes at Xj = L — may not necessarily guarantee 
that its desired value is produced at Xj = L — and t = T — at- However, before studying the value 
of this last plaquette, one needs to fix the value of all plaquettes located at 0 < Xj < L — and 
t = T — at- These plaquettes evaluate to 

T’(o,i)(xi,t = T-at) = U^^^^\xi,T- 

X[/(SQCD)t(x^ + as,T- _ ^^y-ieQAi{^,,T-at)as 

— g-2eQ[Ai(xi,r-at)-Ai(xi,T)]asgjeQ[Ao(xi,T-at)-Ao(xi+as,r-at)]atp(QCD)^^^ j, _ 

X g-ieQ[Ai(xi,r)-Ai{xi,0)]as ^^ 2 '^ 

where we have suppressed the Xj and x^ dependences of the functions. To eliminate the unwanted 
phase factor, one can modify the value of the link along the 0 direction as following 

[/(j^‘^°^(x t) X g*eQAo(x,i)at ^ gieQ[Aj(x,t=0)-Ai(x,i=r)]/o,i(x)x<5gT_,jj 
where the boundary function /o,i(xj) satisfies 

Ai{xi + as,Xj,Xk,t = 0) - Ai{xi + as,Xj,Xk,t = T) /o,i(xi -h as,Xj,Xk) 

= Ai{x,t = 0) - Ai{x,t = T) {fo,i{x) + as). (34) 
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The dependence of the new functions /o,i on any coordinate variable Xj must be realized through 
(xj — [^] L) as before. As will be shown shortly, although this condition on /o,j guarantees that the 
expected values of plaquettes at 0 < Xj < L — and t = T — at are produced, it is not sufficient in 
general to ensure that the desired value of plaquettes is also produced in the 0 — j plane with j ^ i 
when t = T — at- We will return to this point below. 


With the modihcations of the Uq and Ui links according to Eqs. (30) and (33), we are ready to 


inspect the value of the plaquette located at the far corner of the lattice with x* = L — Ug and t = T — at- 
This plaquette evaluates to 

T’( 0 ,i)(T -as,T- at) = - as,T - a^)e*eQAo(L-a.,T-at)atgieQ[A,(L-a., 0 )-Ihi-»s,r)]/o.di^-a.) 

X^{QCD)(^ _ Q'^^ieQAi{L-as,0)as^ieQ[Ao(0,0)-Ao{L,0)]fifi{0) 

X^^QCD)t(o rp _ ieQAo(0,r—at)atg—ieQ[Ai(0,0)—Ai(0,T)]/o_i(0) 

X — 0,^ T — “s,T—at)asg—ieC3[Ao(0,T—at)—Ao(L,T—at)]/i,o('r—at) 

_ g-*eQ[Ai(L-as,T-at)-Ai(L-a3,T)]as+teQ[Ao(T-as,T-at)-Ao(T,T-at)]at x ()x- t) 

, (35) 


X 

L ds ^ ^ 

^ieQ[Ai{xi,0)-Ai{xi,T)\as 

'T—at _ _ ' 

g-*eQ[Ao(0,t)-Ao(T,b]at 


.Xi=0 

. t=o 


where we have used Eqs. (31) and (34) to arrive at the final expression. As is required, the value of the 


plaquette is independent of the fifl and /o,t functions. In order for this plaquette to have the desired 
value, one can impose the condition 


= 1, (36) 


L ds ^ ^ 

TT g-ieQ[Ai{xi,,Xj ,Xfc ,t=0)-Ai{xt,Xj,Xk,t=T)]aa 


'T—at ^ ' 

gieQ[A o (xi=0,Xj,Xk,t)-Ao(xi=L,Xj,Xk ,t )]at 

.Xi=0 


. t=o 


to set the extra phase factors in Eq. (35) to unity. This condition, and its implication for the allowed 


gauge field choices, require further discussion. Let us first point out that Eq. (36) arises from adding 


up the value of the U{1) plaquettes that are corrected by introducing modified links throughout the 
lattice in the 0 — i plane. 


L—CLs T—dt 

n n [< 

x ^=0 t =0 


-ieQ[Ai{xi,Xj,xt;,t)-Ai{xi,Xj,xt.,t+at)]as+ieQ[AQ{xi,Xj,xt.,t)-Ao{xi+as,x.j,Xk,t)]at 


L — ds ^ _ 

TT g-*e<3 l^i (Xi ,x.j ,Xfc ,0) - Ai (Xi ,Xj ,xt.,T)]as 

'T—at _ _ ' 

gieQ[Ao (0,Xj ,XJ, ,t)-Ao{L,Xj .x^,t)] at 

.Xi=0 

. t=o 


_ ieQ [- Exjdo [a (xi ,Xj ,Xk,0)-Ai (xi ,Xj ,Xfe ,T)] as+X)Lo'^* [^O (0,Xj ,Xfe ) - Aq (T,Xj ,Xfe ,4)] at j 


(37) 


This latter exponent is ieQ times the net flux of electric held through the 0 — i plane. Explicitly in 
the continuum limit. 


$ 


rT rL 

|^](“®*(xj,Xfc) = ^ dt J dxj Ei(xi,Xj,Xfc,t) 


dxi 


Aj(xi,xj,xfc,0) - Aj(xj,xj,xfc,r) 


+ / dt 

Jo 


Ao{0,Xj,Xk,t) - Ao{L,Xj,Xk,t) 


(38) 


So the condition in Eq. ( |36[ ) states that the net hux of the electric held through 0 — i plane must 
be quantized, <I>|^|(”'^*(xj,x^) = with n G Z. This condition can only hold if the left-hand side 
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is independent of Xj and coordinates. It therefore constrains the class of electric fields that can 
be implemented in this relatively general framework on a periodic lattice.^ As we will see below, the 
same conclusion can be drawn by inspecting the consistency of the conditions on the functions 
(^, 1/ = 0,1,2,3). 

Let us not constrain our discussion to the classes of gauge fields that can be quantized through 


Eq. (36) for the moment, and continue to allow a general coordinate dependence for the gauge fields. 


We will obtain all such constraints that arise on the space-time dependence of gauge fields shortly. 
Having studied the case of electric field in the i direction in great detail, it is easy to deduce the 
modified links as well as the QC for the scenario where electric field is nonvanishing along other spatial 


directions. Firstly, Eq. (33) must be generalized as following to incorporate additional links in the 
{j 7 ^ i) directions, 

t) X g*eQ^o(x,t)at ^ gieQ[Ai(x,t=0)-Ai(x,i=T)]/o,i(x)X(5t,T-at ^ 39 ^ 

i=l,2,3 


At a first glance, it might appear that a similar equation to that in Eq. (30) must represent the value 


of the link in the i direction with i = 1,2,3. However, one should bear in mind that in this general 
scenario there exist a nonzero magnetic field and the value of the links in the i direction can only be 
fully fixed after inspecting the value of plaquettes in all cartesian planes. 

The value of a plaquette in the i — j plane with 0 < Xj < L — and 0 < Xj < L — is 


_ ieQ[Ai{xi,Xi)-Ai{xi,Xj+as)]aa-ieQ[Aj{xi,Xj)-Aj{xi+aa,x.j)]as 'niQC'D)^ % 


(40) 


which correctly accounts for the magnetic field flux through the surface area aj in the i — j plane in 
the continuum limit. 


4> 


(B) 

(iJ) 


rXi+as 

= --€ijk I dx' 


fXj+aa 


dx'Bfc(x',x'-,Xfc,t) 


'Xi 


= [Ai(xi,xj,xfc,t) - Ai{xi,Xj + as,Xk,t)] 0,3 - [Aj{xi,Xj,Xk,t) - Aj{xi + as,Xj,Xk,t)] as + O (a^) . 

(41) 


Eq. (40) suggests that the results for the modified links, and the QC that was obtained from the 


consideration of the plaquettes in the 0 — i plane, can be simply carried over to the case of the i — j 
plane by substituting 0 ^ i, i ^ j, at ^ ag and T —?■ L in those relations. In particular, one finds that 

jy(QCD)^^^ X g*eQAi{x,0as X g*eQ[Ao(xi=0,Xj,xj,,t)-To(xi=L,Xj,xj,,i)]/i,o(xj,Xfe,t)X(5x^^L_as 

X g*eQ[Aj(xi=0,Xj,Xfc,t)-Aj(xi=L,Xj,Xfc,t)]/ij(xj,xj,,Ox(5xi,i-<i3|j42^ 

where the boundary function fij{xj) must satisfy 

Aj{xi = 0,Xj + as,Xk,t) - Aj{xi = L,Xj + as,Xk,t) + as,Xfc,t) = 

Aj{xi = 0,Xj,Xk,t) - ^j(xj = L,Xj,Xk,t) + as ), (43) 


® More general frameworks that allow more complicated space-time dependences for the fields may potentially exist. 
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Aj{xi = 0,xj,xfc + as,t) - Aj{xi = L,Xj,Xk + as,t) /ij(xj,Xfc + as,t) = 

Aj{xi = 0,Xj,Xk,t) - Aj{xi = L,Xj,Xk,t) (44) 


Aj{xi = 0,Xj,Xk,t + at) - Aj{xi = L,Xj,Xk,t + at) fi,j{xj,Xk,t + at) = 

Aj{xi = 0,xj,xk,t) - Aj{xi = L,Xj,Xk,t) /i,j(xj,Xfc,t). (45) 

While the first condition arises from requiring the plaquettes in the i — j plane to have their desired 
value when xt = L — Ug, the last two conditions arise from setting the value of plaquettes in the i — k 
and i — 0 plane to their desired values. Furthermore, by studying carefully the value of the plaquettes 


in the i — j and i — k planes when Xj = L — and given the modihed links in Eq. (42), one arrives 
at the following conditions on /j^o 

^o(xi = 0,Xj + as,Xk,t) - ^o(xi = L,Xj + as,Xk,t) /i,o(xj + as,Xk,t) 

Ao{xi = 0,Xj,Xk,t) - Ao{xi = L,Xj,Xk,t) /i,o(xj, Xfc, t), (46) 


^o(xi = 0 ,Xj,Xfc + as,t) - ^o(xj = L,Xj,Xk + ag,t) /i,o(xj,Xfc + as,t) 

Ao{xi = 0,Xj,Xk,t) - ^oi^i = L,Xj,Xk,t) /i,o(xj,Xfc,t). (47) 


These conditions add to the condition in Eq. (31) and must be satisfied simultaneously to obtain a 


consistent solution for fifl. Similarly, if the value of the plaquette in the 0 — j and 0 — k planes are 
considered when t = T — at, one arrives at 

Ai{xi,xj + as,Xk,t = 0) - Ai{xi,xj + as,Xk,t = T) fo,i{xi,Xj + as,Xk) 

Ai{x, t = 0) - Ai{x,t = T) /o,i(x), (48) 


Ai{xi,xj,xk + as,t = A) - Ai{xi,xj,xk + aa,t = T) /o,i(xi,x^,x^) 

Ai{x, t = A)-Ai{x,t = T) fo,i{x), (49) 


which supplement the previous condition on /o,j in Eq. (34). 


With these modified links, one can see that, in addition to the three QCs in Eq. (36) for i = 1, 2, 3, 
three more QCs arise 


= 1, (50) 


by requiring that the desired value of the plaquette located at Xj = L — and Xj = L — a* is 


L ds ^ 


L Qis ^ 

TT g-*e(3 [Ai (Xi,Xj =0,Xfc ,t)-Ai (Xi ,Xj =I/,Xfc ,t)]as 


TT [Aj (xi =0,Xj ,Xfc ,t)—Aj {xi=L,Xj ,Xfe ,i)] 

.Xi=0 


J-io 


generated. Eq. (50) is the statement that the net flux of the magnetic held through the i — j plane 
must be quantized, t) = with n G Z. Therefore, if the hux is dependent on the Xk 

{k 7 ^ i,j) and t coordinates, this condition can not be satished in general. 

Equations above for the functions must be satished simultaneously to ensure that the proper 
values of the elementary plaquettes are produced throughout the lattice. However, it is not guaranteed 
that, for any given A held, these equations possess consistent solutions for each To clarify this 
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point, consider the fi^ function which must be obtained recursively from Eqs. (31), (46) and (47). It 
is straightforward to see that these equations are consistent with each other only when depends 
solely on the Xj and t coordinates. In general, there exists a valid only if Ai, solely depends on 
Xf^ and Xi, coordinates. Note that if Ai, is independent of the coordinate, no discontinuity occurs 
in the value of plaquette in the fj, — u plane when = T and = at for fi = 0 while 

for /X = i). As a result no needs to be introduced to guarantee periodicity.® 


Lfj_ = L and = a. 


Interestingly, all such conditions on the space-time dependence of can be shown to be consistent 
with the statement that the net electric or magnetic flux through any plane on the four-dimensional 
lattice (a closed surface in the toroidal geometry) must be space-time independent. This is exactly the 
condition we deduced by examining the QCs in Eq. (36) and (50). In the next section, we present 
the setup for the implementation of several chosen background fields on a periodic lattice and and will 
specify the corresponding QCs. 


IV. EXAMPLES: PERIODIC IMPLEMENTATION OF SELECTED BACKGROUND 
FIELDS ON A HYPERCUBIC LATTICE DEDUCED FROM THE GENERAL CASE 


The examples that follow provide a setup for the implementation of selected background electric and/or 
magnetic fields that preserves the periodicity of the lattice calculation. These are deduced from the 
general case of the previous section, the results of which will be summarized here for convenience. In 
order for the background 17(1) gauge links to be implemented periodically, they must be introduced as 




[ iLu — an 


(51) 


where /x and v assume distinct values. If each depends only on x^ and x^ coordinates, there exist 
functions that satisfy the following recursive relation on the lattice. 


Ai/{x^ — 0, Xp ap'j Ap{x^ Xp -1- ap^ f T ^u) — 

Api^Xfj, = 0, Xp') Ap{xfi = LtLi Xp') {ffi^pi^Xp) Up). (52) 


is vanishing if Ap is independent of x^J Under the conditions specified, the net electric or magnetic 
flux through any plane in the continuum limit is constant and must be quantized. On the lattice these 
QCs read 


1 

P 


Ljy> - CLl/ ^ 

TT ^-ieQ[Afj,{xp,Xp=0)-Ap{Xfj,,Xp=Lp)\ap 


^ieQ[Ap{xp=0,Xp)-Ap(xp=Lp,Xp)^ap 

1- 

II \ 

o 

1_ 


. Xp=0 


(53) 


® Another way to ensure that the desired values of the plaquettes are produced adjacent to the boundaries of the lattice 
is to enforce a set of micro QCs. These QCs can be deduced by setting the extra factor in the value of plaquettes 
near the boundaries equal to one (without requiring any gauge link to b e modified). For example, one can set the 

eQ[Ao(L,Xj,Xfc,t)-A(](O.Xj,Xfc,t)]at | 29 | and 


coordinate-dependent factors e" 


,-.Xfe,T)-Ai(xi,Xj,Xfc.O)]<i, 

Eq. (32 I equal to 1, such that not only the correct values of plaquettes in the 0 —i plane at 0 < < L —Os, t = T —at 

and Xi = L — Oa, 0 < t <T — at are produced, but also the correct value of plaquette at the far corner of the lattice, 
i.e., at Xi — L — as, t = T — at, is produced, without the emergence of an extra factor as in Eq. ( |35[ ). Clearly, the extra 
factor in Eq. | |35[ ), which is a product of all the coordinate-dependent phase factors above, is automatically equal to 1 
because of the micro QCs. Unfortunately, this procedure will not alway work due to the coordinate-dependence of the 
new conditions (note that even if the gauge fields are chosen to be independent of the transverse coordinate, the QCs 
still carry a longitudinal coordinate dependence). On the lattice, where space and time coordinates are discretized, 
gauge fields with simple rational dependences on the coordinates could allow such micro QCs to be satisfied. However, 
such conditions are more restrictive on the background field parameters than the QC derived in this section on the 
total flux of the field, leading to large quanta of background fields that is not desired for most applications. 
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All the space-time dependences in the following examples must be understood through the floor- 
function prescription of Sec. [TT| 


Example I: A constant electric field in the X3 direction 


We have already discussed this case in Sec. [TT] Here we choose a different gauge than that taken in 
Sec. 0 


= (Ao, -A) = (0,0,0, ait ), 


(54) 


which produces an electric field in the X3 direction, E = 01X3. The nontrivial gauge links that produce 
this background field are 


jjiQCD) jj{QCD) ^ g—ieQair/o,3(x3)x<5t^T-at 

^(QCD) ^(QCD) ^ ^ieQaitas 


(55) 

(56) 


where /o ,3 satisfies /o, 3 (x 3 TOs) = /o, 3 (x 3 ) TOs, with the solution /o, 3 (x 3 ) = X 3 once one sets /o, 3 ( 0 ) = 
0.® There is only one QC in this case, 


L ds 

^-ieQaiTds 

X3=0 


^—leQaiTas 


as 


_ ^—ieQaiTL 


1 , 


which constrains the value of the field strength, oi, 

27rn 

oi = —^-, 

eQTL 

with n G Z. 


(57) 


(58) 


Example II: A constant magnetic field in the X3 direction 


One can pick the following gauge potential 

Af, = (Ao, -A) = 1 ^ 0 , yX2, -yxi, 0 ^ , 


(59) 


to generate a uniform magnetic field in the X 3 direction, B = 61 X 3 . The nontrivial gauge links that 
are required to implement this background field are 


f/WCD) ^|QCD) ^ g|eQfeiX2as ^ ^^eQt'iLfi^2{^2)xSxi,L-as 
jjiQCD) ^^QCD) ^ ^—ieQfeixias ^ g—fe(36iL/2,i(xi)X(5x2,L-as 


(60) 

(61) 


^ In cases where does not depend on the x^, coordinate, as is the case in most of the examples in this section, the 
satisfies the relation f^^^ixi, + = fij.,vixv) + o-a with the solution = x^, once one sets f,j.,u{0) — 0. If 

Ai, depends on both and Xi, coordinates, it is possible to transform to a gauge where Ai, does not depend on x,^, 
as long as the condition = 0 is satisfied, where is the EM field strength tensor and no summation over 

/i and V is assumed. 

® Note that the initial value of the function is arbitrary as it drops out of the value of plaquettes. 






18 


The conditions on /i _2 and / 2 ,i are respectively, /i, 2 (x 2 + a^) = /i, 2 (x 2 ) + a* and / 2 ,i(xi + o^) = 
/ 2 ,i(xi) + tts, with the solutions /i ^2 = ^2 and f 2 ^i = xi, once the initial values of the functions are set 
to zero. There is one QC for this choice of the field 


L CLs 

gfeQbiLas 

.xi=0 


L CLs 



gfeQbiLas 

.X2=0 


^ieQbiLds 


_ ^ieQbiL^ _ ^ 


which constrains the strength of the magnetic held parameter, 


(62) 


with n G Z. 


bi 


27rn 

egL2’ 


(63) 


Example III: A space-dependent magnetic field and a constant electric field 

When interested in extracting the spin polarizabilities of nucleons, one can choose the following 
gauge potential 


= (- 4 o, -A) = (0, 62 x 1 x 2 ,0, a2t ), 


(64) 


to produce a space-dependent magnetic field, B = 62 X 1 X 3 , as well as a constant electric field, E = 02 X 3 . 
These background fields generate a nonvanishing interaction proportional to 20 'i(ViBj+VjBj)Ej (with 
CTj denoting Pauli matrices) in the nonrelativistic Hamiltonian of the spin-| hadron in external fields, 
which gives access to the 'yEiM 2 spin polarizability of the hadron, see Ref. [22]. To produce this electric 
field on a periodic lattice, the nontrivial gauge links to be implemented are 


^(QCD) ^ g—ieQa2r/o,3(x3)X(5t_'r-ai 


rr(QCD) 

^(QCD) ^ jj{QCD) ^ gieQb2XiX2as 


a 

u. 


(QCD) 

2 

(QCD) 


'1 

^^QCD) ^ g-ieQ62XiL/2,i(xi)X(5x2,L-as 
^(QCD) ^ ^ieQa2tas 


(65) 

( 66 ) 

(67) 

( 68 ) 


The /o ,3 and / 2 ,i functions satisfy recursive relations /o, 3 (x 3 + as) = /o, 3 (x 3 ) + a* and / 2 ,i(xi + as) = 
(/2,i(xi) + a^), respectively, with solutions /o, 3 (x 3 ) = X 3 and /2,i(xi) = for xi > 0 , once 

the initial values of the functions are set to zero. The only QCs are 


L ds 

^ieQb2XiLas 

xi=0 


gfeQ 62 D"(D-ad ^ gfeQb2D3(l-^) ^ 


(69) 


and therefore 


L CLs 

^-ieQa2Tas 

X3=0 


^—ieQa2Tas 


L_ 

as 


— g—ieQci2T L 


1 , 


62 


47rn 

eQL3(l-^)’ 


02 


2Trn' 

eQTL' 


(70) 


( 71 ) 


with n, n! G Z. Note that if we had only quantized the flux of the magnetic field in the continuum 
limit, we would have introduced a deviation from periodicity on the lattice oi 0{^). To avoid such 
discretization errors one must quantize the fields according to Eq. (531, thereby respecting the lattice 
geometry. 
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Example IV: A time-dependent electric field 


Another spin polarizability of nucleon can be accessed via the background gauge potential 


= (^ 0 , —A) 


0 , -ast , 04^,0 


( 72 ) 


which produces a time-dependent electric field, E = 03 ^x 4 + 04 X 2 .® This background field generates 
a nonvanishing interaction proportional to cr • E x E in the nonrelativistic Hamiltonian of the spin-^ 
hadron in external fields, which gives access to the 'JEiEi spin polarizability of the hadron, see Ref. 
|22) . To produce this electric field on a periodic lattice, the nontrivial gauge links to be implemented 
are 


jjiQCD) jj{QCD) ^ g-|eQa3r2/o,i(xi)x<5t^T-ot ^ g-*eQa4T/o,2(x2)X(5t,T-at 

^^(QCD) ^(QCD) ^ ^^eQast^as 

jj(QCD) _^ ^(QCD) ^ ^ieQa^tas 


(73) 

(74) 

(75) 


The /o,i and /o ,2 functions satisfy recursive relations /op(xi + a^) = /o,i(xi) + and fo, 2{^2 + Os) = 
/o, 2 (x 2 )+ a<j, respectively, with trivial solutions /o,i(xi) = xi and /o, 2 (x 2 ) = X 2 , once the initial values 
of the functions are set to zero. The only QCs are 


n 




xi=0 


gleQasT^as 


gieQasT^L _ ^ 


(76) 


and therefore 


with n, n' G Z. 


Z/ Qjs 

^ieQa4Tas 

X 2=0 


^jgQqj/^T Qjs 


as 


_ ^ieQa4TL 


= 1, 


dvrn 2Trn' 

^ 5 - 7 5 

eQT‘^L eQTL 


(77) 


(78) 


Example V: A space-dependent electric field and a constant magnetic field 

A third spin polarizability of nucleon can be accessed with the following gauge potential 

\ = (Ao, -A) = (^-yx^, 0,0, 63 x 1 ) , (79) 

which produces a space-dependent electric field, E = 05 X 2 X 2 and a constant magnetic field, B = 63 X 2 . 
This background field generates a nonvanishing interaction proportional to 2 < 7 i(ViEj■+VjEi)Bj in the 
nonrelativistic Hamiltonian of the spin-| hadron in external fields, and so gives access to the 'JM 1 E 2 


® Any time-dependent magnetic field will necessarily generate a time-dependent flux through the transverse plane. This 
means that one can not implement a periodic background fleld with the method of this work to access the spin 

polarizability of hadron. This is, up to a numerical factor, the coefficient of the cr ■ B x B term in the nonrelativistic 
Hamiltonian of the spin-| hadron in an external magnetic field. 
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spin polarizability of the hadron, see Ref. j22j . To produce this electric field on a periodic lattice, the 
nontrivial gauge links to be implemented are 


r(QCD) 

0 

. rr(QCD) 
Uq 

X 

-|eQa5x|at 

(80) 

dQCD) 

-7 

X 

g-2eQb3L/i,3(x3)X(5x^,L-as 

(81) 

r(QCD) 

2 

. rr(QCD) 
U 2 

X 

/ 2,0 ^ ^xo — O.S 

^ 1 

(82) 

r(QCD) 

3 

. rr(QCD) 
U 3 

X 

^ieQbz'x.ias 

(83) 


The /i ^3 and / 2 ,o functions satisfy recursive relations /i, 3 (x 3 + a^) = /i^ 3 (x 3 ) + and / 2 ,o(^ + CLt) = 
f 2 ,o{t) + CLt, respectively, with trivial solutions /i, 3 (x 3 ) = X 3 and f2,o{t) = i-, once the initial values of 
the functions are set to zero. The only QCs are 


T—at 

n 


^-^eQasL'^at 


t=0 


1 

^-CeQasL^at “t _ ^-^eQa^L'^T _ ^ 


(84) 


and therefore 


with n, n' G Z. 


L Ojg 

^-ieQbsLas 

X3=0 


^—ieQbsLas 


as 


_ ^—ieQbsL'^ 


1 , 


“5 


dvrn 

eQL^T' 


bs 


27rn' 

eQL^' 


(85) 


( 86 ) 


Example VI: A plane-wave electric field 

As suggested in Ref. |39) . background EM plane waves can be used to calculate the off-forward 
matrix elements of current operators between hadronic states, enabling an extraction of hadronic form 
factors.Additionally, a recent proposal in Ref. m demonstrates the advantage of a plane-wave 
background field in evaluating the hadronic vacuum polarization function with lattice QCD. This 
approach proceeds by constraining the polarization function using the susceptibilities with respect to 
the background magnetic field amplitude at specific momenta. Due to the condition on the space-time 
dependence of the flux of the background field in each plane, our periodic implementation of EM plane 
waves will be limited to fields with certain Fourier modes. For example, an electric field of the form 
E = e®‘^'^X 3 with q* 7 ^ 0 for all i = 1, 2, 3 will generate a coordinate-dependent flux. We can however 
generate an electric field of the form, e.g., E = a 6 e*'^®^®X 3 , from the following gauge potential 

= (Ao,-A) = 0 , 0 , 0 ) , (87) 

with a constant flux through the 0 — 3 plane. Since the form factors in the continuum (infinite-volume) 
limit are rotationally invariant, the constraints on the components of the momentum transfer vector in 
this setup will not prevent one from accessing the form factors at any values. The only limitation on 


In general, any space-time inhomogeneity in the background field gives access to such off-forward quantities due to 
injecting energy and/or momentum. Plane-wave background fields have the advantage of isolating contributions with 
a given momentum transfer. 
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the (magnitude) of the transferred momenta may arise from the implementation of fields on a periodic 
lattice as will be deduced below. 

To produce the electric field chosen above on a periodic lattice, the nontrivial gauge links to be 
implemented are 


U, 


(QCD) 


0 


U. 


(QCD) 




( 88 ) 

(89) 


where /s^o satisfies the recursive relation / 3 ,o(i + flt) = fsflit) +®i with trivial solution / 3 ,o(i) = t, once 
the initial value of the function is set to zero. The only QC is 


T—at A 

e <?3 '' * 

t =0 


e 13 


a-t eQa^ q ^io-^L\rp 

= e 13 ^ ^ = 1. 


For any arbitrary value of the field amplitude parameter, as, this QC can be satisfied with 


qs = 


27rn 

ir^ 


(90) 


(91) 


with n G Z. This constraint on means that the additional U{1) phase factor in Eq. (89) is equal to 


unity. It also means that with a background field method, the EM form factors can only be accessed 
at quantized values of momentum transfer, the situation which is also encountered when form factors 
are calculated through a direct evaluation of hadronic matrix elements on the lattice. However, one 

(r) 

could allow for nonqunatized values by placing conditions on the real part, Og , and imaginary part, 
ag\ of ag. Indeed, by requiring 


(j) _ nq^n' (r) _ sin(g3L) (j) 


dc. — 


eQT 


Ojc - 


1 - cosiq^L) 


(92) 


with n' £ and q^ ^ for n G Z, the QC in Eq. (90) is satisfied. As we already saw, for q^ = 
the QC in Eq. (90) trivially holds.The result of this latter case is indeed consistent with the periodic 


implementation of an oscillatory magnetic field (through sin and cos functions) in Ref. |l0] in which 
the Fourier modes of the applied field are taken to be quantized. 


V. CONCLUSION 

The introduction of classical electromagnetic fields in lattice quantum chromodynamics calculations 
provides a powerful technique to obtain a variety of electromagnetic properties of hadrons and nuclei. 
To extend the utility of this technique, the current implementations of uniform background fields 


A plane-wave background gauge field as in Eq. (871 has a nonvanishing imaginary piece. It therefore results in a non- 
Hermitian fermionic determinant (after integrating out the quark fields) which can hinder the probabilistic evaluation 
of the associated path integral (analogous to the effect of a nonvanishing chemical potential). For isovector quantities, 
one can avoid this sign problem imposed by such plane-wave gauge field by only implementing the background fields 
on the valence quarks. For isoscalar quantities, where there are disconnected contributions to the matrix elements 
of the EM current, contributions from a charged quark sea can not be ignored. For these quantities, where a full 
imposition of the background fields on the valence and sea quarks is required, the computational cost associated with 
a non-Hermitian measure in the Monte Carlo sampling of the path integral can be controlled by tuning the amplitude 
of the gauge field to be small. Note that the amplitude of the field can be made arbitrarily small (but nonvanishing) 
only if qs is quantized according to Eq. (911. Alternatively, one can implement only real oscillatory functions, i.e. 
sinlgaXs) or cosiqaxa), to access the desired off-forward matrix elements. 
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|22[ 1231123H3S] have been extended in this paper to the case of background fields that have temporal 
and/or spatial nonuniformities. Such field configurations provide access to static and quasi-static 
properties of hadrons and nuclei, such as their higher EM moments and polarizabilities as well as their 
charge radii. They also provide means to directly extract EM from factors as energy and momentum 
are injected into the hadronic system immersed in appropriate external fields |39) . Such possibilities can 
be explored in upcoming lattice QCD calculations once the corresponding background fields required 
for each quantity is consistently implemented on the particular lattice geometry. Periodic boundary 
conditions are widely used in lattice QCD calculations, and given that such boundary conditions 
typically result in a simpler hadronic theory in general to be matched to lattice QCD calculations, 
it is important to perform the background field calculation with the imposition of these boundary 
conditions. In this work, we have considered the most general space-time dependent U{1) gauge fields 
imposed on QCD gauge configurations, and have shown that under certain conditions on the space-time 
dependence of the fields, a periodic implementation of background U{\) gauge fields is possible. 

To make the discussions more transparent before getting into formalities of the general case, we 
have first presented the special case of an electric field generated through a scalar potential with an 
arbitrary dependence on one spatial coordinate. The necessary modifications to the naive U{1) links 
adjacent to the boundary of the lattice are obtained by ensuring the expected values of the elementary 
plaquettes are produced throughout the lattice. These expected values are nothing but the 17(1) 
phases corresponding the flux of the electric field in the continuum limit through surface areas at x as 
at each space-time point on the lattice, where at and a* refer to lattice spacings in temporal and spatial 
directions, respectively. Additionally, a quantization condition is obtained that ensures that the flux 
of the electric field is quantized. From this special case, two examples of a uniform, and a linearly 
varying electric field in space, are constructed. We have numerically confirmed that only the periodic 
implementation of gauge fields, according to the prescription proposed, gives rise to smooth correlation 
functions for neutral pions across the boundary of the lattice. 

For the general case of gauge fields with arbitrary space-time dependences, one can follow the 
same procedure as that of the special case above, except that obtaining the modified links near the 
boundary is more involved. This is simply because of the fact that when a component of the gauge 
field depends on more than one space-time coordinate, fixing the values of plaquettes in one plane to 
their desired values (by modifying the links adjacent to the boundary) can potentially affect the values 
of the plaquettes in other planes. By carefully accounting for such possibilities, we have derived a 
set of equations for the functions that need to be introduced in the modified links and have discussed 
the conditions on classical fields that guarantee the existence of solutions for these equations. We 
have further shown that these conditions are equivalent to the statement that the flux of the electric 
and magnetic fields through each plane of the lattice must be coordinate independent and quantized. 
The latter is a condition that must be met to ensure the expected value of the plaquette at the far 
corner of the lattice in each plane is produced. In a parallel approach, we have shown with details 
in the appendix that these conditions arise from a more general class of boundary conditions, namely 
electro/magneto-periodic boundary conditions |131 HHHTTl 1501453] . where one assumes that the gauge 
fields are only periodic up to a gauge transformation. 

We have used our general construction to explicitly work out several examples relevant to the 
extraction of various EM moments, spin polarizabilities and form factors from future lattice QCD 
calculations. Within our construction, time-dependent magnetic fields can not be studied in this 
framework, which limits access to the 'JMiMi spin polarizability of the nucleon. However, it is plausible 
that a more general construction will allow for this polarizability to be extracted. A rather interesting 
case is a plane-wave background EM field, which can be devised to be periodic on a lattice by properly 
choosing the Fourier modes of the fields. Through studying the corresponding QC, it is apparent that 
to require periodicity, either the Fourier modes in each direction must be quantized for arbitrary values 
of the field amplitude parameter (as implemented in Ref. m), or the Fourier modes can be chosen 
to be arbitrary while the amplitude parameter must be quantized. Our results have applications in 
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upcoming lattice QCD calculations that aim to extract such quantities using the background field 
method. 
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APPENDIX: ELECTRO/MAGNETO-PERIODIC BOUNDARY CONDITIONS AND THE 
ASSOCIATED CONDITIONS ON THE BACKGROUND FIELDS 

The QCs obtained in this paper for the parameters of the background fields, and the conditions that al¬ 
lowed a periodic implementation of C/(l) gauge links on the lattice, can also be deduced from imposing 
more general boundary conditions, namely electro/magneto-PBCs. These boundary conditions require 
the gauge and matter fields to be periodic up to a gauge transformation so that all gauge-invariant 
quantities will be periodic. The electro/magneto-PBCs have been introduced, and extensively dis¬ 
cussed, by’t Hooft in Refs. US SS EO] for the case of Abelian and non-Abelian gauge theories, and 
were later adopted to derive the QC for the case of uniform background fields implemented on a torus 
Ha HUMl EU Eg. Here, we aim to make explicit the relation between’t Hooft’s conditions and those 
presented in this work.^^ In particular, we obtain those electro/magneto-PBCs that give rise to the 
same conditions on the U{1) gauge fields that have been obtained in this paper with PBCs. It is shown 
that the conditions that are placed on the gauge functions when electro/magneto-PBCs are imposed 
(see below) are equivalent to the conditions on the flux of the EM field when PBCs are imposed. 
The discussions presented in this section correspond to the continuum spacetime and so can only be 
compared with the continuum limit of the QCs presented in Eq. ( |53[ ) (i.e., the conditions on the flux 
of the field strength tensor in the continuum limit). With a lattice geometry, one needs to use the 
results presented in earlier sections of this paper. 

Consider a gauge field that depends on all space-time coordinates, Xjy, Xp, Xo-), with all 

indices being distinct.For the moment, let us assume that only the component of the gauge 
field is nonzero. We demand that the field is periodic up to a gauge transformation at the boundary. 
Explicitly, 


Xi/^ Xp, X(j) 

— ^pi Xcr) ~\~ (Xp,, X^, Xp, Xfj), 

(93) 

Xp, Xfj) 

— -Ap(xp, 0, Xp, Xfj) ~\~ ^(xp, Xp, X(j); 

(94) 

Xy, Lp^ Xfj) 

— Ap(Xp, Xy^ 0, Xcr) ~\~ (Xp, Xj/, Xfj), 

(95) 

Xyi Xp, L(j) 

= Ap(xp,x^,Xp,0) + 9pf7(^’'^)(xp,Xi.,Xp). 

(96) 


In order for the D functions to represent a gauge transformation, they must transform the matter field 
Tp at the boundaries as well. Since the transformation of the matter fields depends on H (and not only 


For another approach in elucidating the connection between these boundary conditions for the case of a uniform 
magnetic field on the lattice, see Ref. |48j . 

The coordinate-dependence of the functions in this appendix must not be realized through the floor-function prescrip¬ 
tion. In fact, all the functions are tilde functions as defined in the main text, but the tilde over functions is dropped 
to keep the notion cleaner. 
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the derivate of Q with respect to x^), relations (93][%) do not entirely fix the boundary conditions on 
the matter fields. We now show that the following choice of electro/magneto-PBCs on matter fields, 


■lp{Lp,Xy,Xp,Xa) 

= 'ljj{0,Xy,Xp,Xa), 

(97) 

P{Xp,Ly,Xp,X„) 


(98) 

1p{Xp,Xy,Lp,X„) 


(99) 

1p{Xp,Xy,Xp,L„) 


(100) 


along with the boundary conditions on the fields above, gives rise to the same QCs on the back¬ 
ground fields’ flux as obtained in this paper by imposing PBCs. First note that the transformations 
in Eqs. (93 96) and (97||100) are consistent with the following solutions for the gauge functions. 


(np rp rp T* ^ 

= Xp [Ap{Lp, Xy, Xp, Xa) - 

■ Ap{{),Xy,Xp,Xa)] , 

(101) 

Q^^’’''>{xp,Xp,Xa) = 

/ dx' [Ap{Xp,Ly,Xp,Xa 
Jo 

) - Al^{x'p,Q,Xp,X^)] , 

(102) 

Ll^^'^\Xp,Xy,X„) = 

/ dx'[Ap{x',Xy,Lp,x„ 
Jo 

) - Ai^{Xp,Xy,Q,Xa)] , 

(103) 

II 

fXp 

/ dx'[Ap{Xp,Xy,Xp,L^ 

Jo 

) - Al^{x'p,Xy,Xp,Q)] . 

(104) 


Now consider the transformation of field ^|J at the corner of the lattice in the n — v plane, i.e., at 
and Xy = Ly. One may first transform the '0 field with function and then with the 

function. 


Ly, Xp, Xa) = 

or the other way around, 

Ly, Xp, x^) = Q, Xp, x^). 

These two relations are compatible if 


(105) 


(106) 


(107) 


which can only hold in general if is independent of the Xp and x^ coordinates. From Eqs. (93) 

and (94), these conditions should hold if Ap_ is independent of these coordinates. Now given that from 
Eq. (102), n(^’’^)(0, Xp, X(j) = 0 and Xp, Xq-) is the total flux of the EM field through the 

fj, — v plane on the lattice (recalling the assumption that only Ap is nonvanishing), the condition in 


Eq. (107) is identical to the continuum QC as obtained in Sec. Ill 

rLu. 


'0 


dx p\Ap(x p, Ly') Mp(xp,0)] — 


27rn 


n G Z. 


(108) 


It might seem that by alternatively considering the gauge transformed field at the corner of the 
lattice in any other plane, one could arrive at the same condition which would be eliminating any 
space-time dependence of the gauge field. However, it is easy to see that this is not the case. For 
example, by considering the gauge transformed matter field in the corner of the lattice in the p — a 
plane one arrives at 


iplx 




X, 


M ^p1 ^ 


,0)^ieQn(i^’P^Xp,x 


„L, 




/, 0 , 0 ), 


(109) 
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as well as, 




but it is not hard to see from Eqs. (103) and (104) that 

+ 9!^^"'P\xp,Xy,L„) = 9!^^"'P\xp,Xy,Q) + QS^’''\xp,,Xu,Lp). 


( 110 ) 


( 111 ) 


As a result, the two transformations for ip{xp,Xu, Lp, L„) are identical and no extra conditions must 
be placed on the Xp and x^ dependences of the Ap field. 

Now let us assume that a second component of the gauge field, A^y, is also nonvanishing. Then 
the same equations as those in (93]|%) and (101 104) can be written for the gauge transformed A^, 
field at the boundaries, and for the corresponding solutions for the gauge functions, respectively. The 
transformations of the field turn into 


'tp{Lp,Xu,Xp,Xa) 

= Xu, Xp, Xa] 

), 

(112) 

'lp{Xp,Lv,Xp,Xa) 

= 0, Xp, Xy, 

), 

(113) 

'lp{Xp,Xiy,Lp,Xa) 

_ {Xfj. ,Xv ,Xa) (xu,Xu, 

''"‘^^'il;{xp,Xu,0,Xo-), 

(114) 

'lp{Xp,Xiy,Xp,Lyj) 

_ ^ieQVl£’'^'>{xu,Xu,Xp) ^ieQQ.^t^’'^'>(xu,Xu 

’^p^'il;{xp,Xu,Xp,0), 

(115) 


with our choice of boundary conditions on the matter field. The same argument discussed above based 
on the compatibility of the gauge transformations of the matter fields at the corner of the lattice in 
each plane now leads us to the condition 


{0,Xp,Xa)—U(f^’'''> {Lfj,,Xp,Xa)] ieQ[U^''’i^'> (0,Xp,Xa)—i^^''’>^\Liy,Xp,Xay^ _ 


(116) 


which requires that both the Ap and Ai, fields to be independent of the Xp and x^ coordinates. 
Additionally, it is easy to see that this condition implies that the total flux of the EM field through 
the fj, — u plane is quantized, 

27m 

/ dxp[Ap{xp,L^)-Ap{xp,0)]- dx^[A^{Lp,x,y) - Aiy{0,Xu)] =n £ 7., (117) 

Jo Jo eQ 

conditions that were all deduced in Sec. |III| by imposing PBCs. In general, if any of the Ap component 
of the gauge field with /r = 0,1,2,3 is nonvanishing, the compatibility relations similar to the ones 
considered above must be studied at the corner of the lattice in all planes, and they are all satisfied if 
the flux of the field strength tensor through each plane is coordinate-independent and quantized. 

We conclude this appendix by a remark. Note that the boundary conditions in this section are 
defined by only identifying the point Lp with the point 0 in the /r direction, and similarly in other 
directions. If one imposes more restrictive conditions such that all points Xp -|- Lp are identified with 
points Xp, the QCs derived in this section for the flux of the background field through the ^ — u plane 
will depend, in general, on the Xp and x^ coordinates and can not be quantized. This is in contrast with 
the case of uniform background fields where both the PBCs and electro/magneto-PBCs can be satisfied 
with a single space/time-independent QC. Such system possesses a translational invariance in units of 
Lp for all /X (corresponding to the magnetic translational group in the case of a uniform magnetic field, 
see Ref. m)- With a background EM field that is space-time dependent, no translational invariance 
exists prior to imposing the boundary conditions. Therefore, the lack of discrete translational invariance 
in units of Lp should not come as a surprise. Given our choice of periodifying the functions with the 
use of the floor function, we have however explicitly built up such translational invariance in the setup 
presented in this paper. One may however wonder whether this choice can be distinguished from a 
choice for which only the point Lp is identified with the point 0 on the lattice. This answer is that with 
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a lattice action that at most depends on the first derivative of the gauge and matter fields, these two 
choices are identical at the practical level. Explicitly, one only evaluates fields at points 0 < 
and only the values of the fields at must be specified through the boundary conditions. 
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